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)FOURI FT ft TRANSFORM PkOGKAM IN PDP-b FORTRAN II 
i NTNO IS THE NUMBER OF FUNCTION POINTS 
IFNTNO IS NTNO IN FLOATING POINT 
IN IS NT JO/2 (FRACTIONAL PART IS TRUNCATED) 

IN INDEXES THE HARMONICS 
IIX INDEXES THE FUNCTION VALUES 

IK (=IX +1) INDEXES THE STORAGE OF FUNCTION VALUES 
IFX IS IX IN FLOATING POINT 
IA IS ThE COSINE TERN COEFFICIENT 
10 IS THE SINE TERN COEFFICIENT 

ilj (=S*FX*T) IS THE ARGUMENT OF THE SINE AND COSINE 
IFUNCTIONS. IT IS COMPUTED IN THREE STEPS TO MINIMIZE 
I UNNECESSARY COMPUTATIONS 


15,10 


VALUES 


I READ NUMBER OF 
I LESS THAN 201 
I TYPE 105 
I ACCEPT 110,NT WO 
I IF (NTNO-200) 15, 

I TYPE 115 
i GO TO 5 
i N=NTN0/2 
;FNTNO=NTNO 
;DI MENSI ON F(200) 

I READ IN FUNCTION 
IDO 20 K=1,NTNO 
;IX=K-1 
iTYPE 120,IX 
;ACCEPT 125,F(K) 
t CONTINUE 
JTYPE 130 

JS=2.0*3.141593/F NTNO 
;ITERATE OVER HARMONIC 
»DO 35 M=0,N 
;A=B=0.0 
I FM = M 
; T=FM*S 
l J=NTWO-1 

I SUM OVER FUNCTION VALUES 
i DO 30 IX=0,J 
1FX=IX 
i U= T*FX 
i K=IX + 1 

;A=A+F(K)*COSF(U) 

} B=B+F < K)*SINF(U) 

I CONTINUE / cf.<. «. 

JIF CM) 305,32,305 
i IF (NTN0-2+M) 31,32,31 
IA=2•0+A/FNTWO 

;b=2.0*r/fntwo 

;TYPE 135,M,A,B 
)GO TO 35 
)A=A/FNT»0 
iB=0.0 

l TYPE 135,M,A,B 
)CONTINUE 
l PAUSE 


FUNCTION POINTS AND CHECK THAT NUMBER 
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IH1T CONTINUE KEY TO 

RESTART P k 0' ■ 


IGO TO 5 


105 

i FORMATC /, "NUMBER OF 

DATA POINTS 

110 

I FORMAT( I ) 


115 

IFQKMATC'NUMBER TOO 

LARGE") 

120 

S FORMAT<"DATA POINT 

",I," = ") 

125 

iFORMATC E) 


130 

i FORMAT!/,"HARMONIC 

COSINE 

135 

JFORMAT!I," ",E,E,/) 


i end 
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